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WEIGHTED SPANNING TREES ON SOME SELF-SIMILAR GRAPHS 

DANIELE D'ANGELI AND ALFREDO DONNO 

Abstract. We compute the complexity of two infinite families of finite graphs: the Sierpiiiski graphs, 
which are finite approximations of the well-known Sierpihsky gasket, and the Schreier graphs of the 
Hanoi Towers group H^ 3 ' acting on the rooted ternary tree. For both of them, we study the weighted 
generating functions of the spanning trees, associated with several natural labellings of the edge sets. 

^H ! 

O 
(N 

£h ' 1- Introduction 
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The enumeration of spanning trees in a finite graph is largely studied in the literature, and it has many 
applications in several areas of Mathematics as Algebra, Combinatorics, Probability and of Theoretical 
Computer Science. 

Given a connected finite graph Y = (V(Y), E(Y)), where V(Y) and E(Y) denote the vertex set and 
the edge set of Y , respectively, a spanning tree of Y is a subgraph of Y which is a tree and whose vertex 
r) | set coincides with V(Y). 

The number of spanning trees of a graph Y is called the complexity of Y and is denoted by t(Y). The 
famous Kirchhoff's Matrix- Tree Theorem (1847) states that t(Y) is equal to (the constant value of) any 
cofactor of the Laplace matrix of Y, which is obtained as the difference between the degree matrix of Y 
and its adjacency matrix. Equivalcntly, t(Y) ■ \V(Y)\ is given by the product of all nonzero eigenvalues 
of the Laplace matrix of Y. 

It is interesting to study complexity when the system grows. More precisely, given a sequence {T„}„>i 
of finite graphs with complexity r{Y n ) 1 such that |V(Y^)| — > oo, the limit 

hm l ^m, 

O : |v(y„)|-K» \V(Y n )\ 

when it exists, is called the asymptotic growth constant of the spanning trees of {Y n } n >i (see |13j). 
[~^. ■ A spanning k-forest of Y is a subgraph of Y which is a /c-forest, i.e., it is a forest with k connected 

^^ ' components, and its vertex set coincides with V(Y). 

The enumeration of spanning subgraphs, in general, for a graph Y, is also strictly related to the Tuttc 
polynomial Ty{x,y) of the graph: more precisely, it is known that Ty (1,1) equals the complexity of 
Y, Ty(2,l) equals the number of spanning forests of Y, and Ty (1,2) is the number of its connected 
spanning subgraphs (see [5J [§] , where this analysis is developed for the finite Sicrpihski graphs and for 
other examples of finite graphs associated with the action of automorphisms groups of rooted regular 
cH | trees). 

A finer invariant of the graph Y is a finite abelian group <&(Y), whose order is exactly the complexity 
of Y. This group occurs in the literature under different names, depending on the context. It was 
introduced in [T] as the Picard group of Y (or the Jacobian of Y), whereas it is shown in [3] that the 
Picard group is isomorphic to the group of critical configurations of the chip-firing game on Y. As any 
finite abelian group, $(Y) can be decomposed into direct sum of invariant factors. The dependence of 
this decomposition on the properties of Y has been studied by several authors, (see, e.g., [12]), but not 
much is known so far. Explicit computations have been performed for certain families of graphs. 

In many optimization problems it is often useful to find a minimal spanning tree of a weighted graph. 
Hence, it is interesting to study spanning trees when a weight function on E(Y) is introduced. So let 
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w : E(Y) — > R + be a weight function defined on the edges of Y. Let T be the set of all spanning trees 
of Y. With each spanning tree t G T, we can associate the weight 

W(t):= J] w(e), 

eG-E(t) 

i.e., the product of the weights of the edges of Y belonging to E(t). Then, the weighted generating 
function of the spanning trees of Y is the polynomial on the variables w(e), e G E(Y), given by 

T(w):=Y,W(t). 
ter 
It follows from the definition that, when evaluated on the constant weight w = 1, the generating function 
yields the complexity of the graph, since in this case one has W(t) — 1, for each t G T. 

In this paper, we will study weighted spanning trees on two infinite families of finite graphs very close 
to each other: the Sicrpihski graphs, which are finite approximations of the famous Sierpinski gasket, 
and the Schrcicr graphs of the Hanoi Towers group H 1 - 3 ' , which is an example of a self-similar group (see 
Definition 11.21 below) . 

We recall some basic facts about self-similar groups. Let T q be the infinite regular rooted tree of degree 
q, i.e., the rooted tree in which each vertex has q children. Each vertex of the n-th level of the tree can 
be regarded as a word of length n in the alphabet X = {0, 1, . . . , q— 1}. Now let G < Aut(T q ) be a group 
acting on T q by automorphisms generated by a finite symmetric set of generators S. Suppose, moreover, 
that the action is transitive on each level of the tree. 

Definition 1.1. The n-th Schreier graph S„ of the action of G on T q , with respect to the generating 
set S , is a graph whose vertex set coincides with the set of vertices of the n-th level of the tree, and two 
vertices u, v are adjacent if and only if there exists s G S such that s(u) = v. If this is the case, the edge 
joining u and v is labelled by s. 

The vertices of £„ are labelled by words of length n in X and the edges are labelled by elements of S. 
The Schreier graph is thus a regular graph of degree |5| with q n vertices, and it is connected since the 
action of G is level-transitive. 

Definition 1.2 ( |14j). A finitely generated group G < Aut(T q ) is self-similar if, for all g G G,x G X , 
there exist h G G,y G X such that 

g(xw) = yh(w), 
for all finite words w in the alphabet X. 

Self-similarity implies that G can be embedded into the wreath product Sym(q) I G = Sym(q) x G q , 
where Sym(q) denotes the symmetric group on q elements, so that any automorphism g G G can be 
represented as 

9 = a(go,...,g q -i), 
where a G Sym(q) describes the action of g on the first level of T q and gi G G,i — 0, ...,q — 1, is the 
restriction of g on the full subtree of T q rooted at the vertex i of the first level of T q (observe that any 
such subtree is isomorphic to T q ). Hence, if x G X and w is a finite word in X , we have 

g(xw) =a{x)g x {w). 

The class of self-similar groups contains many interesting examples of groups which have exotic prop- 
erties: among them, we mention the first Grigorchuk group, which yields the simplest solution of the 
Burnside problem (an infinite, finitely generated torsion group) and the first example of a group of in- 
termediate growth (see [TU] for a detailed account and further references). In the last decades, the study 
of automorphisms groups of rooted trees has been largely investigated: R. Grigorchuk and a number 
of coauthors have developed a new exciting direction of research focusing on finitely generated groups 
acting by automorphisms on rooted trees [3] . They proved that these groups have deep connections with 
the theory of profinite groups and with complex dynamics. In particular, for many examples of groups 
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belonging to this class, the property of self-similarity is reflected on fractalness of some limit objects 
associated with them [13]. 

Since the Schreier graphs are determined by group actions, their edges are naturally labelled by the 
generators of the acting group and it takes sense to study weighted spanning trees on them, with respect 
to this labelling. 

The paper is structured as follows. In Section [2j we study weighted spanning trees on finite approxi- 
mations of the well-known Sicrpihski gasket, endowed with three different edge labcllings: 

• the "rotational- invariant" labelling, whose special symmetry allows to explicitly compute the gen- 
crating function of the spanning trees (Theorem 12 -2[) and to perform a statistical analysis about 
the number of edges, with a fixed label, occurring in a random spanning tree of the graph (Propo- 
sition gTTJl; 

• the "directional" labelling, where the weights depend on the direction of the edges; for this model, 
the weighted generating function of the spanning trees is described via the iteration of a polyno- 
mial map ( Theorem 12. 8|) : 

• the "Schreier" labelling, strictly related to the labelling of the Schreier graphs of the Hanoi Towers 
group H^ 3 >; also in this case, the weighted generating function of the spanning trees is described 
via the iteration of a polynomial map (Theorem I2.12p . 

In all these models we follow a combinatorial approach. The self-similar structure of the graph (in 
the sense of |16j ) allows to study both unweighted and weighted subgraphs recursively More precisely, 
we introduce three different generating functions associated with spanning trees, 2-spanning forests, 3- 
spanning forests and, using self-similarity, we are able to establish recursive relations (Theorems 12.11 12.61 
and I2.10p and to give an explicit description of them (Theorems 12.21 12.81 and I2.12|) . More generally, the 
self-similar structure of a graph turns out to be a powerful tool for investigating many combinatorial and 
statistical models on it: see, for instance, [71 151 1151 [TrJl [TT] . 

In Section [31 we consider the Schreier graphs of the Hanoi Towers group H^\ whose action on the 
ternary tree models the famous Hanoi Towers game on three pegs (see |llj). endowed with the natural 
edge labelling coming from the action of its generators. Even if these graphs also have a self-similar 
structure, the combinatorial approach used in the case of the Sierpihski graphs seems to be much harder 
here. Therefore, our technique consists in using a weighted version of the Kirchhoff's Theorem: we 
construct the Laplace matrix by using the self-similar presentation of the generators of the group, which 
is impossible in the case of Sierpihski graphs, where there is no group structure. In this case, the 
generating function is described in terms of iterations of a rational map (Theorem 13. 5|) : this kind of 
approach already appears in [21 [UJ (see also [5] , where we use the same strategy to compute the partition 
function of the dimer model on the Schreier graphs of the Hanoi Towers group) . 

2. Spanning trees on the Sierpinski graphs 

The problem of enumeration of spanning trees in Sierpihski graphs was largely treated in literature 
(see, for instance, [SI [H]). We consider here three different labcllings of the edges of these graphs and 
write down the associated generating function of the spanning trees. In all the models, the self-similarity 
of the graphs plays a crucial role to study the problem recursively. The description of the generating 
function strongly depends on the symmetry of the labelling of the graph: as we will sec, in the first model 
that we consider, which is invariant under rotation, we are able to give an explicit formula for it; in the 
two remaining models, where we do not have invariance under the action of any symmetry group, the 
generating function is described via the iteration of two polynomial maps. 
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2.1. First model: "Rotational-invarianf'labelling. Let F\ be the graph in the following picture. 




For each n > 1, we define, by recurrence, the graph T„ + i as the graph obtained by partitioning an 
equilateral triangle in four smaller equilateral triangles and by putting in each corner a copy of T n . 
Observe that this labelling of the graph is invariant with respect to the rotation of ^ . We represent in 
the following picture the graph T 2 . 




We want to study weighted spanning trees on the graphs {r„}„>i. For each n > 1, we put: 

• T n (a, b, c) = weighted generating function of the spanning trees of r„; 

• S n (a 7 b, c) = weighted generating function of the spanning 2-forests of r„, where two fixed outmost 
vertices belong to the same connected component and the third outmost vertex belongs to the 
second connected component; 

• Qn(a, b, c) — weighted generating function of the spanning 3-forests of r n , where the three out- 
most vertices belong to three different connected components. 

Observe that, because of the rotational invariancc of the labelling of the graph, the function S n (a,b,c) 
does not depend on the choice of the two outmost vertices. In what follows, we will often omit the 
argument (a, b, c) of the weighted generating functions. 

Theorem 2.1. For each n > 1, the weighted generating functions T n {a,b,c),S n {a,b 1 c) and Q„(a, 6, c) 
satisfy the following equations: 



(1) 


T n +i — 6T n S n 


(2) 


"n+l = <T n o n + T n Q n 


(3) 


Q n+1 = 12T„S n Q n + USl 


with initial conditions 





Ti(a, b, c) = 3(a + b)(ab + ac + be) 2 
Si(a,b,c) = (a + b)(a + b + 3c) (ab + ac + be) Qi(a,b,c) = (a + b)(a + b + 3c) 2 . 

Proof. The graph I\ l+ i can be represented as a triangle containing three copies of r„. 
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Wc will use the pictures 






to denote, respectively, the case where in a copy of r„: 

• the three outmost vertices are in the same connected component; 

• two outmost vertices are in the same connected component and the third one is in a different 
connected component; 

• the outmost vertices are in three different connected components. 

The only way to construct a spanning tree of r„ + i is to choose a spanning tree in two copies of T n and 
a spanning 2-forcst in the third one, as in the following picture. 




This argument proves Equation (TTJ), where the factor 6 is given by symmetry (we have to take into 
account both reflections and rotations). 

Next, we are going to prove Equation @ (we analyze, for instance, the case where the leftmost and 
the rightmost vertices are in the same connected component). Consider the two following pictures. 





These possibilities, together with their symmetric, obtained by reflecting with respect to the vertical axis, 
give a contribution to SVi+i equal to 4T n S^. Consider now the two following configurations. 





Since the picture on the left has to be considered together with its symmetric, we get a contribution to 
5,1+1 equal to 3T n S^. Finally, the contribution T^Q n is described by the following picture. 




D. D'ANGELI AND A. DONNO 



This proves Equation ©. 
We have now to prove Equation (|3]) about Q n +i- Consider the following situations. 





They provide, by symmetry, a contribution equal to 12T n S n Q n - The following pictures give, by symmetry, 
a contribution of 125^ to Q n +i- 





Finally, the following two pictures give a contribution of 25^ to Q n +i- 





This completes the proof. 



□ 



Theorem 2.2. For each n > 1, the weighted generating functions T n {a,b,c) 1 S n {a,b 1 c) and Q n (a,b,c) 
satisfying Equations fTj). @ and ([3]), with the initial conditions given in Theorem \2.1[ are: 



T n (a, 6, c) = 2" 



-1 3^ + 271-1 3 re-i_ 2Tt+1 Q „_ i , 3^- + -1 , s 3 TI + 1 

_ 3 3 5 i (a + fr^ (a + b + 3c) = (a& + ac + bc)^~ 



S n (a,b,c) =2" 5 3 J 5 r~ (a + 6) 3 (a + 6 + 3c)^ (a& + ac + be) ~^~ ; 
Q„(a, 6,c) = 2 s 3 a 5 3 (a + 6) 3 (a + b + 3c) s (a& + ac + bc)^~ . 

Proof. The proof is by induction on n. It is easy to verify that, for n = I, one gets the initial conditions 
given in Theorem 12.11 Then, one can check that the functions given in the claim satisfy Equations dXJ) , 
([2]) and (|3|). We omit here the explicit computations. □ 



It follows that T„(l, 1, 1) = r(r n ); similarly, s(T„) := S n (l, 1, 1) is the number of spanning 2-forests of 
r„, where two fixed outmost vertices belong to the same connected component and the third outmost 
vertex belongs to the second connected component; q(T n ) := Q n (l, 1, 1) is the number of spanning 
3- forests of T n , where the three outmost vertices belong to three different connected components. 

Corollary 2.3. For each n > 1, one has: 

, , 3"-l 3"+ 1 +2re+l 3 n -2n-l 

(1) r(T„) = 2—3 * 5^^; 

, 3"-l 3 T1 + 1 -2n-3 3"+2ra-l 

(2) s(r„) = 2—3 — 3 — 5^^; 

, , , , 3"-l 3" + 1 -6ti-3 3"+6ti-1 

(3) q{T n ) = 2~ 3 * 5^^. 

In particular, the asymptotic growth constant of the spanning trees of T n is 3 log 2+2 log 3 + g log 5 . 
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Proof. It suffices to evaluate the weighted generating functions described in Theorem l2.2l for a = b = c = 
1. The asymptotic growth constant is then obtained as the limit 



lim 



log(r(r n )) 



where \V(T n )\ = §(3 n 



\V(T n )\ ' 
1) is the number of vertices of r„, for each n > 1. 



□ 



Remark 2.4. The same values of the complexity and of the asymptotic growth constant have been found 
in [6] and JT5] , where the authors study unweighted spanning trees of T n . 



2.2. Second model: "directional" labelling. Consider now a new sequence of graphs {r„}„>i, which 
coincide, as unweighted graphs, with the graphs studied in Section 12.11 and whose edges are endowed 
with a new labelling, that we call directional labelling. It is clear that an edge of r„ can point in three 
different directions: up (from left to right), down (from left to right) or horizontal. Then, we label by 
a each edge pointing up, by b each horizontal edge, and by c each edge pointing down, where, as usual, 
a,b,c € K.+ . Here we draw the three first examples. 




Remark 2.5. Observe that the indices are now shifted by 1 with respect to the case of the rotational- 
invariant labelling considered in Section [2. II the reason is that a rotational-invariant labelling using three 
labels a, b and c cannot be defined on a simple triangle. 



In this section, we study the weighted spanning trees of the graph T r 
labelling. For each n > 1, we put: 



endowed with the directional 



• T n (a, b, c) = weighted generating function of the spanning trees of r„; 

• U n (a, b, c) = weighted generating function of the spanning 2-forests of r„, where the leftmost and 
the rightmost vertices belong to the same connected component, and the upmost vertex belongs 
to the second connected component. Similarly, by rotation, we define R n (a,b,c) (respectively 
L n (a,b,c)) for the spanning 2-forests of T n , where the rightmost (respectively leftmost) vertex is 
not in the same connected component containing the two other outmost vertices; 

• Qn{a, b, c) = weighted generating function of the spanning 3-forests of T n , where the three out- 
most vertices belong to three different connected components. 

Observe that, in this model, we need to introduce three different functions U n (a 7 b,c),R n (a,b 7 c) and 
L n (a,b,c), since the edge labelling is not invariant with respect to a rotation of ^ as in the previous 
case. On the other hand it is clear that, for each n > 1, one has U n (l, 1, 1) = i? n (l, 1, 1) = L n (l, 1, 1) 
and this common value is equal to S n _i(l,l,l), where S n (a,b,c) is the generating function introduced 
in Section |2"7T1 In what follows, we will often omit the argument (a, b, c) of the generating functions. 
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Theorem 2.6. For each n > 1, the weighted generating functions T n (a,b,c), U n (a,b,c), R n (a,b,c), 
L n (a,b,c) and Q n (a,b,c) satisfy the following equations: 

(4) T n+1 =2T 2 (U n + R n + L n ) 

(5) U n+1 = T n U n (2R n + 2L n + 3U n ) + T 2 Q n 

(6) R n +i = T n R n (2L n + 2U n + 3R n ) + T n Q n 

(7) L n+1 = T n L n (2R n + 2U n + 3L n ) + T 2 n Q n 

(8) Qn+i = 4T n Q n (U n + Rn + L n ) 

+ 2 (U 2 {R n + L n ) + Rl(L n + U n ) + Ll(R n + [/„)) 
+ 2U n R n L n , 
with initial conditions 

Ti(a, 6,c) = ab + ac + be U\{a,b 1 c) = b i?i(a, 6, c) = a Li(a, 6, c) = c Qi(a, 6, c) = l. 

Proof. It is easy to check that the initial conditions hold. Then, the proof of each recursive equation 
follows the same strategy as in Thcorcm l2.ll □ 

Remark 2.7. Observe that, by replacing U n , R n and L„ with S n , one finds again the equations given 
for the rotational-invariant model in Theorem 12.11 



In order to get explicit solutions of the equations given in Theorem 12.61 we put 

<pi (a, b, c) = ab + ac + be <p2 (a, 6, c) = a + b + c 
/(a, 6, c) = 3a~b + 3ab" + 3a c + 3ac + 36 c + 36c + 7a6c 
and let us define the function F : M 3 — >1 3 as F(x,y,z) = [F\{x,y,z),F2{x,y, z),F^(x,y,z)), where 
F\{x,y,z) = 3x 2 +3xz+3xy+yz F2(x,y,z) = 3y 2 +3xy+3yz+xz Fs(x,y,z) = 3z 2 +3xz+3yz+xy. 

Moreover, we denote by F± (a, b, c) the i-th coordinate of the vector F^ = F(. . . F(F(a, 6, c))), where 
the function F is iterated k times. Note that F^ '(a, 6, c) = Fi(a, 6, c), for each i = 1, 2, 3. 
Finally, for each k > 3, put <^>fc(a, 6,c) = <f)^_\{F\{a, b, c),F2(a, b, c), F^{a, b, c)), so that 

fe ( a ,6, C ) = 2 (F(' £ - 2 )(a,6,c)). 

Theorem 2.8. The weighted generating functions T n (a, 6, c), U n (a,b,c), R n (a,b,c), L„(a, 6, c) and 
Q n {a, 6, c) satisfying Equations ((4]), ([5]), ([6]), ([7]) anrf ([8]), wii/i the initial conditions given in Theorem, 



l +.-j 



T„(a, 6, c) = 2 12" II <£fc 6 (<*,&, c )j for each n> I; 

fc=i 

n — 1 n _ fe m 

U n (a,b,c) = 2 it 1^ 8 (a,6, c)F2 Il_1) (a, 6, c), for each n> 2; 

fc=i 

i?„(a,6,c) = 2 3 i=" +3 11^ 5 (a,6,c)F 1 ( "' 1) (a,6,c), /or eadn > 2; 

fc=i 

n— 1 n _(.ij 

L„(a, 6, c) = 2 3 "~i'" +3 II f k B (a, 6, c)F 3 (ti_1) (a, 6, c), /or eac/i n > 2; 



k=\ 
n-1 



t (a,6,c) ^2 3 " 1 i2" +39 H^ fc 6 (a,6,c)/(F (Tl - 2) (a,6,c)) for each n> 3, 



fe=i 



wii/i C/i(a, 6, c) = 6, i?i(a, 6, c) = a, Li(a, 6, c) = c, Qi(a, 6, c) = 1 and Q2{&, 6, c) = 2/(a, 6, c). 
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Proof. The proof works by induction on n. One can directly find: 

U 2 {a,b,c) — 0i (a, b, c)F 2 (a, b, c) Ri{a,b,c) = <f>\(a, b, c)i<i(a, b, c) L 2 {a,b,c) = </>i(a, b,c)F 3 (a,b,c) 

Ti(a,6,c) = 4>i(a,b,c) Q 3 (a,b,c) = 24>l(a,b,c)f(F 1 (a,b,c),F 2 (a,b,c),F 3 (a,b,c)), 

and so the basis of induction holds. We only prove the assertion for T n (a, b, c), by showing that Equation 
((4]) is satisfied (the computations in the other cases are similar but more complicated). One has: 



_+5 _ 3 n -6ra+3 

(o,6,c)-2— n— 



2T%(U n + R n + L n ) = 2 • 2 — ^^ J| ^ fc 3 (a, b, c) • 2^^ _[]_ </) fc 6 (a,6,c) 

fc=i fe=i 

■ (F 2 ( "- 1) (a,6,c)+F 1 "- 1) (a,6, C ) + ^"- 1) (a,6, C )) 



3 Ti-r-i +6T1 _3 



n 



fc=i 



(a, 6, c) 
(jf- 1} (a, 6, c) + F[ n ~ l) (a, b, c) + F^ (a, b, c)) 

" + 1 .". 3 ?i-fc + 2 +3 

2 3 """ 3 J] <t>~ 6 (a, b, c)4>2 [F^-V (a, b, c)) 
fe=i 



3 71 ^ 1 +6n-3 



ra+1 

fc=i 



- -I- ^ 



(a,6,c) = T )l+ i 



a 



2.3. Third model: the "Schreier" labelling. Consider the graph I?i in the picture below and define 
by recurrence, for each n > 1, the graph T n+ i as constituted by the union of three copies of T n in the 
following way: for each one of the outmost vertices of T n+ i, the corresponding copy is given by the graph 
r„, reflected with respect to the bisectrix of the corresponding angle. 




Remark 2.9. For each n > 1, the graph T n coincides with the graph obtained from the Schreier graph 
S n of the Hanoi Towers Group H^ by deleting the three loops and by contracting all the edges joining- 
two different elementary triangles, keeping the labels in £„ on the remaining edges. (See Section [3.11 ) 

Define the generating functions T n (a, b, c), U n (a, 6, c), R n (a, b, c), L n (a, b, c) and Q n (a, b, c) to have the 
same meaning as in the case of the directional labelling (Section l2.2p . In what follows, we will often omit 
the argument (a, 6, c) of the generating functions. 
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Theorem 2.10. For each n > 1, the weighted generating functions T n (a,b,c), U n (a,b,c), R n (a,b,c), 
L n (a,b,c) and Q n (a,b,c) satisfy the following equations: 

(9) T n+1 = 2Tl (U n + R n + L n ) 

(10) U n+1 = T n (3L n R n + U n R n + U n L n + 2U%) + T*Q n 

(11) R n +i = T„ (3U n L„ + U n R n + R n L n + 2R n ) + T n Q n 

(12) L„+i = T„ (3U n R n + L n U„ + R n L n + 2L 2 n ) + T%Q n 

(13) Qn+i = iT n Q n (U n + Rn + L n ) 

+ 2 (U 2 {L n + R n ) + R 2 n {U n + L n ) + Ll(R n + U n )) 
+ 2U n R n L n , 
with initial conditions 

Ti(a,b,c) = ab + ac + be Ui(a,b,c) = b R\{a,b 1 c)=a L\{a,b,c) — c Qi(a, o, c) = 1. 

Proof. It is easy to check that the initial conditions hold. Then, the proof of each recursive equation 
follows the same strategy as in Theorem 12. II □ 

Remark 2.11. Observe that, by replacing U n , R n and L n with S n , one finds again the equations obtained 
for the rotational- invariant model in Theorem 12.11 

Put 

tpi (a, b, c) = ab + ac + be tj)2{(l,b,c) = a + b + c 

g(a, 6, c) = 3a~b + Sab + 3a c + 3ac + 3b c + 36c + labc 
and let us define the function G : M 3 — > R 3 as G(x,y,z) = (Gi(x,y,z),G2{x,y,z),Gs(x,y,z)), where 
Gi(x, y,z) = x +2yz + xy + xz G2(x,y, z) = y +2xz + xy + yz Gs(x,y,z) = z +2xy + xz + yz. 
Finally, for each k > 3, put ipk(a, b, c) = V'fe-i(Gi(a, 6, c),G2(a, b,c), Gs(a, b,c)), so that 

Ma,b,c)^^j 2 (G^- 2) (a,b, C )y 

Theorem 2.12. The weighted generating functions T n (a,b,c), U n (a,b,c), R n (a,b,c), L n (a,b,c) and 
Q n (a,b,c) satisfying Equations §§§, ([TO"]) , (jTTJ) , (|12p and (H2J), with the initial conditions given in Theorem 
WM are: 

T n (a,b,c) = 2 2 I I i() k 2 (a,b,c) for each n > 1; 
fe=i 

)i-i 



3 n 1 l -m w- 3 1 / -. \ 

U n (a,b,c) = 2 a 11 ^fe 2 ( a ^? c )^2 ~ ( a A c ) / or eac ^ n > 2 ; 
fc=i 

n— 1 rt— fe 

O Tl 1 -^ _ _ 3 — 1 / -I \ 

i?„(a,o,c) = 2 3 II V^ 2 (a, 6, c)G 1 "~ (a, 6,c) for each n > 2; 
fc=i 

n — 1 n_ fc 

L n (a, 6,c) = 2 s I I V'fe 2 ( a )fr) c )G*3 (a, 6, c) /or eac/i n > 2; 



fc=i 

n-2 



Q„(a,6,c) = 2 £L ^- L J] Vfc 2 3 (a,6,c).grG ( "" 2) (a,o,c)") foreachn>3, 



fc=i 



wii/i t/i(a, 6, c) = 6, i?i(a, 6,c) = a, Li(a, o, c) = c, Q/i(a, 6, c) = 1 and Q2(a, b, c) = 2g(a, 6, c). 
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Proof. The proof is by induction on n. One can directly find: 

U 2 (a,b,c) = 2-0i (a, b, c)G 2 (a, b, c) R 2 (a,b,c) = 2^i(o, b, c)Gi(a, b, c) L 2 (a,b,c) = 2ipi (a, 6, c)G 3 (a, 6, c) 

Ti(a, 6, c) = V>i(a,6,c) Q 3 (o,6,c) = 2 4 ipl(a, b, c)g(G x (a, b,c),G 2 (a, b,c),G 3 (a, b, c)), 
and so the basis on the induction holds. We only prove the assertion for T n (a, b, c), by showing that 
Equation (J5|) is satisfied (the computations in the other cases are similar but more complicated). One 
has: 

n n-l_ U ~ 1 3"- fc -l 

2T%(U n + R n + L n ) = 2-2 3 "" " _i n^r +1 (a,6,c)-2^^^ J] V fc ~^ (a, 6, c) 

fc=i fc=i 

• (G^"- x) (a, 6, c) + Gi n-1) (a, 6, c) + G^ (a, 6, c)) 
= 2 ^ n^r~^^ 1 ( a ' 6 ' c )^(G (n - 1) (a,6,c)) 



fe=i 



™+ ] on-fc + 1 



2i ^ 1 II^k 2 KM=T n+1 . 



fc=i 



a 



Remark 2.13. Note that the function g(a,b,c) coincides with the function f(a,b,c), introduced in 
Section 12.21 However, the functions G(x, y, z) and F(x, y, z) do not coincide, which implies that the 
functions ipi(a, b, c) and </>i(a, b, c), factorizing the weighted generating functions of the spanning trees in 
the directional model and in the Schreier model, are different. 

3. Spanning trees on the Schreier graphs of the Hanoi Towers group 

In this section, we study spanning trees on the Schreier graphs of the Hanoi Towers group H^- 3 \ which 
are very similar to the Sicrpihski graphs studied in the previous sections. We start with a combinatorial 
recursive approach as in Section [2J but this leads us to some equations that we can explicitly solve only in 
the unweighted case. Hence, in order to compute the weighted generating function, we follow a different 
approach: we use the self-similar presentation of the generators of the group to write the adjacency 
matrix of the graphs. Therefore, we arc able to write the associated Laplace matrix and then we get the 
weighted generating function of the spanning trees by using a weighted version of Kirchhoff's Matrix- Tree 
Theorem. The self-similarity of the group is the key property that we use. 

3.1. The Schreier graphs of the Hanoi Towers group. Recall that the Hanoi Towers group H^ is 
generated by the automorphisms of the ternary rooted tree having the following self-similar form: 

a = (01) (id, id, a) b = (02) (id, b, id) c= (12) (c, id, id), 

where (01), (02) and (12) are elements of the symmetric group Sym(3), acting on the set X = {0, 1, 2}. 
Observe that a,b,c are involutions. The associated Schreier graphs are self-similar in the sense of [16] . 
that is, £ n _|_i contains three copies of E n glued together by three edges, that we call "exceptional". By 
Definition 11.11 each vertex of £„ corresponds to a word of length n in the alphabet X. The graphs 
{£„}„>! can be recursively constructed via the following substitutional rules |llj . 

11m 

2 M b W lu 
=>■ a/ \c 

20m/ \02m 

b/\ °/\ 

00u a 10ub I2u c 22m 
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22u 



Rule II 




12u« 
b/ 



K02u 
\c 



lOu, 

a/ 



,01u 

VI 



Rule III 



(hi 



Ov 



OOu 



OOw 



Rule IV 



lu 



lv 



llu 



OOu b 20u a 

2 

< 

Rule V a 



21m 



2u 



Uv 



2v 



llu 
22u 



22t> 



Notice that the word u in Rules I and II can also be empty and the words u and v in Rules III, IV, V 

can also satisfy u = v (in this case we get the three loops of £„). The starting point of this recursive 

construction is the Schreier graph Si of the first level. We also draw a picture of £2. 

b 





For each n > 1, the graph £ n +i can be also obtained in the following recursive way: we take the union 
of three copies of £ n and, for each one of the outmost vertices of £„+i , the corresponding copy is reflected 
with respect to the bisectrix of the corresponding angle. 

Remark 3.1. Observe that, for each n > 1, the graph £„ has three loops, centered at the vertices 
0™,1™ and 2", labelled by c, b and a, respectively. This is an easy consequence of the definition of the 
generators a, b and c of H^ 3 \ Moreover, these are the only loops in £„. However, in what follows, we will 
be studying spanning trees on the Schreier graphs £„ without loops, since a spanning tree of £„ cannot 
contain any loop. By abuse of notation, we still denote by £„ the graph without loops. 

3.2. Computation of the complexity. For each n > 1, let T n (a,b,c), U n (a,b,c), R n (a,b,c), L n (a,b,c) 
and Q„(a, b, c) be the generating functions having the same meaning as in Section [2] In what follows, we 
will often omit the argument (a, b, c) in the generating functions. 

Theorem 3.2. For each n > 1, the weighted generating functions T n (a,b,c), U n (a,b,c), R n (a,b,c), 
1 (a, b, c) satisfy the following equations: 

T„4-i = T^{ab + ac+ be) + 2abcT*(U n + R n + L n ) 



L n (a, b, c) and 
(14) 

(15) 



l -n+l 



u, 



71+ 1 — 



ac - 

^3 1 t~i2 



(16) 
(17) 



Rn-\ 



1 — 



J-* n 



+ 1 



bT£ + T,i {(ab + ac+ bc)U n + 2b(aR n + cL n )) 
abcT n {3R n L n + U n (L n + R n + 2U n )) + abcT*Q r , 

aT* + Tl((ab + ac+ bc)R„ + 2a(bU n + cL n )) 
abcT n {3U n L n + R n (L n + U n + 2R n )) + abcT%Q n 

cTl + Tl{{ab + ac + bc)L n + 2c(aR n + bU n )) 
abcT n (3R n U n + L n (U n + R n + 2L n )) + abcT*Q n 



WEIGHTED SPANNING TREES ON SOME SELF-SIMILAR GRAPHS 

(18) Q„+i = 4abcT n Q n (U n + R n + L n ) 

+ T 2 ((2b + a + c)U n + (2a + b + c)R n + (2c + a + b)L n ) 

+ TlQ n {ab + ac + bc) + Tl 

+ 2abc (U 2 {R n + L n ) + R 2 n {U n + L n ) + L 2 n {U n + R n ) + U n R n L n ) 

+ 2T n (U n R n (ac + bc+ 2ab) + U n L n (ab + ac + 2bc) + R n L n (ab + bc+ 2ac) 

+ bU 2 {a + c)+ aRl(b + c) + cL 2 n {a + b)) , 

with initial conditions 

T\(a,b,c) = ab + ac + be U\{a,b,c) = b R\(a,b,c)=a L\(a,b,c) = c Q\{a,b,c) = l. 

Proof. The graph E n+ i can be represented as 



l:s 



S„+i 



where the triangles 1, 2 and 3 represent subgraphs of E„ + i isomorphic to the graph E„. The edges 
labelled by a, b and c in the picture above are the exceptional edges joining the different copies of E„. In 
the pictures of this proof we will use the same conventions as in Theorem 12.11 

Now, a spanning tree of E„ + i can be obtained by choosing a spanning tree for each one of the triangles 
1, 2 and 3, and omitting one of the edges a, b and c in order to have no cycle. This gives the contribution 
Tn{ab + ac + be) to T n+ \. A spanning tree of T, n+ \ can also be obtained by choosing a spanning tree for 
two of the triangles 1, 2 and 3, and a 2-forest in the third triangle. This time, we do not need to omit 
any one of edges a, b and c (see the picture below). 




This situation corresponds to the contribution 2abcT 2 (U n + R n + L n ) to T„+i and Equation (|14[) is proven. 



We want to prove Equation (|15[) . The two following pictures show that a spanning 2-forest of type 
U n +i can be obtained starting from three spanning trees of level n (in this case we omit two exceptional 
edges) , but also by taking two spanning trees in two copies of S„ and a spanning 2-forest of type U„ in 
the third one (by omitting one of the three exceptional edges). 
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More precisely, the first picture corresponds to a contribution equal to bT%, the second one to the 
contribution T%U n (ab + ac + be) to U n+ \. Consider now the following pictures. 





These configurations, together with their symmetric ones obtained by reflecting with respect to the 
vertical axis, give a contribution to [/„+i equal to 2bT%(aR n + cL n ). Consider now the two following 
pictures. 





The left picture, together with its symmetric, gives the contribution 2abcT n R n L n ; the right one, together 
with its symmetric, gives the contribution abcT n U n (R n + L n ). Consider now the two following situations. 





The picture on the left, together with its symmetric, gives a contribution equal to 2abcT n U% to f/ n +i- 
The picture on the right contributes by the summand abcT n R n L n . Finally, the contribution abcT^Qn is 
described by the following picture. 




So we have proven Equation (|15D about U n +\. Equations (|16[) and (|17[) can be proven in a similar way. 
We want to prove now Equation (TTgj) about Q n +i- Consider the following pictures. 





They provide, by symmetry, a contribution equal to AabcT n Q n (U n + R n + L n ). Next, consider the 
following cases. 
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The picture on the left gives, by symmetry, a contribution equal to 2T^(bU n + aR n + cL n ). The picture 
on the right gives T%((a + c)U n + (b + c)R n + (a + b)L n ). Now consider the following cases. 





The picture on the left gives the contribution T% to Q n +i, the picture on the right contributes by 



T*Q n {ab- 



bc). Consider the following pictures. 





By symmetry, each one of the pictures above gives to Q n +i a contribution of 
abc (Ul{R n + L n ) + Rl(U n + L n ) + L 2 n (U n + R n )) . 
Consider now the following pictures. 





They give a contribution equal to 2abcU n R n L n to Q n +\. Now look at the following configurations. 




AiA 




By symmetry, they correspond to the contribution 2T n (abU n R n + bcU n L n + acR n L n ). The following 
picture gives the contribution T n (bU^(a + c) + aR^(b + c) + cL 2 n {a + 6)). 




Hi 
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Finally, we have to consider the following four situations. 





The first and the third pictures provide each a contribution equal to T n (U n R n (a + b)c + U„L n (b + c)a + 
R n L„(a + c)b). The second one gives T„ (bU„(a + c) + aR^(b + c) + ch\{a + b)). The fourth one gives 
2T n (abU n R n + acR n L n + bcU n L n ). This completes the proof. □ 

Equations (fT4")l . (|T5j) , (|TB|) . (fTT)) and (|T51) . with the initial conditions given in Theorem 13. 2[ seem to 
be very hard to be explicitly solved. For this reason, in the next section we use a different strategy to 
find the weighted generating function. On the other hand, if we are interested in the complexity of these 
graphs, we can evaluate each generating function ma = b = c=l and we get simpler equations. Let us 
define: 

• r„ := t(£„) = T„(l, 1, 1) = complexity of £„; 

• s n := U n (l, 1, 1) = -R„(l,l,l) = L n (l, 1, 1) = number of spanning 2-forests, where two fixed 
outmost vertices are in the same connected component and the third one lies in a different 
component; 

• Qn '■— Qn(ljl]l) = number of spanning 3-forests, where each component contains exactly one 
outmost vertex. 

Corollary 3.3. For each n>l, the values T n ,s n and q n satisfy the following relations: 

t 11+ i = 3t^ + 6t^s„ 

Sn+l = Tn + 7 T l s n + 7t„S^ + T^q„ 



q n+ i = 3r n q n + I2r n s n q n + 14s,- 



12r^s„ +rf + 36t„sI 



with initial conditions 



Ti = 3 si = qi = 1. 



Proof. The proof follows from Theorem 13.21 by evaluating each function for a = b = c = 1 and recalling 
that U n (l, 1, 1) = i?„(l, 1, 1) = L„(l, 1, 1), for each n > 1. □ 

Proposition 3.4. For every n > 1, t/ie values r n ,s n and q n satisfying the relations given in Corollary 

I are: 



(1) Tn 

(2) S„ 

(3) qn 






5" -3" 



. (£ ^y 



/n particular, the asymptotic growth constant of the spanning trees of S n is j (log 3 + log 5). 
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Proof. The proof can be given by induction on n. Then, the asymptotic growth constant is obtained as 
the limit 

,. log(r„) 

hm - — 

n->oo |K(2J n )| ' 

where |F(E„)| = 3" is the number of vertices of E„, for each n > 1. □ 

3.3. Computation of the weighted generating function. In this section, wc compute the weighted 
generating function of the spanning trees on the graph E n by using the weighted version of the Kirchhoff 's 
Theorem. The idea is to use the self-similar presentation of the generators of the group in order to describe 
recursively the adjacency matrix of the graph. The Laplace matrix of E ra is then obtained as the difference 
between the degree matrix of E„ and its adjacency matrix. By using the Schur complement Formula, 
one can compute a cofactor of the Laplace matrix passing from a square matrix of size 3™ to a square 
matrix of size 3 n , which turns out to have the same structure of the original matrix, where the entries 
have been transformed via a rational function. This gives rise to a recursion process allowing to compute 
the cofactor. The same strategy is used in [5] to compute the partition function of the dimer model. See 
also [TT], where the authors use the same idea to study the spectrum of the group H^K 

Let A„ be the adjacency matrix of S„, whose rows and columns are indexed by the vertices of E„, i.e., 
words of length n in the alphabet {0, 1,2}, which are ordered lcxicografically. Hence, A„ = a n + b n + c n , 
where the matrices a n , b n and c n describe the action of the generators a, b and c of H 1 - 3 ' , respectively, on 
the n-th level of the rooted ternary tree. In other words, we have: 

J a if a(Vi) = Vj 

\ a n)ij — \ _ , 

I otherwise, 

where Vi, for i = 1, . . . , 3™, denotes the i-th vertex of E„ with respect to the lexicografic order (similarly 
for b n and c n ). Notice that, in this way, we are writing the adjacency matrix of the graph with loops. 
However, in the Laplace matrix, the entries corresponding to loops will be deleted by using the degree 
matrix (a + b + c)I n , according with the fact that a spanning tree cannot contain any loop. 

The matrices a n , b n and c„ can be recursively represented, by using the self-similar description of the 
generators a, b and c of the group. For n = 1, we set 

'c 0^ 

0c 
,0 c 0y 




Cn-l 














cln-l 





cln-l 






Therefore, Kirchhoff 's Theorem states that the weighted generating function of the spanning trees on E r 
can be obtained by computing any cofactor of the Laplace matrix 

((a + b + c)I„_i - Cn-i -al n -i -bl n -i 

-al n -i (a + b + c)I n -i - &n-i -dn-i 

-bl n -i -cJ n -i (a + b + c)I n _i - a n -i, 

We can choose, for instance, to compute the cofactor associated with the first row and the first column 
of the Laplace matrix. In order to compute it we put, for every n > 1, 



Cn-l 


al n-l 


0i n-l 


<-l 


b n -l 


dn-i 


blLi 


dn-1 


Cln-l 
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with 



F - I 



(1 ••• 0\ 

••• 

: '•• 

\0 0/ 



Moreover, we fix the notation jl^ := 7/° + (a + b + c)(Ik — 1°), for each k > 1. 
Define 

((a + 6 + c)I n _i - Cn-l -O^n-l _ ^n-l 

-oJ°_! (a + 6 + c)I„_i - 6„_i -ci„_i 

-W"-i -c/„-i (a + 6 + c)I n _i - a„_i, 

The introduction of the matrices 7° guarantees that e a ;, b is the generating function of the spanning 
trees of £„, because we have performed all the necessary cancellations in A„. More precisely, since 
(A„)n = a + b and this is the only non-zero entry of the first row and column of A„, it turns out that 

c \b IS ec l ua l to the cofactor of (a + b + c)I n — A„ associated with the first row and the first column. 

Moreover, we define the rational function P : R 9 — > R 9 as P(x) = (Pi(x), . . . , Pg(x)), where 

Pl(x)=Xi, P2{x)=X 2l P3(x)=X 3 , 

Pi(x) = l/D 

{x2Xzx\x^XqXq + X\X4x\x 7 x\ — X\x\x\x% — X\x\x^x\x'j + X2X3X4XI + X1X4X5X6X7X8 

— x\x^x^x§ — x\x-2,xzx\ + x\x^xqx^x\ — xix^x^x^Xg — xixIx^xqXjXs — x\x\x^x%xi; 
+ xix 2 x 3 x^xe + X2XzX4X 5 x 6 — xix 3 X4X e x$ + xiX4X 6 x 7 X8 — xix 4 x 6 X7 + X2X3X 4 x^xeXg) 

P 5 (x) = l/D 

(x\Xzx\x\ — X1X2X3X5 — 2:2X4X5X6 + X2X4X5X6X7X9 + X2X4X5X7X9 — X2X4X5X9 
+ X1X3X4X5X6X8 — x\x2x\x$Xj + X2X5X5X7X9 — X2X3X5XgXg + X1X3X4X5X6X8 — X2X5X6X7 
+ X2X4X6X 7 X 9 — X2X 3 X4X6X g — X2X 4X^X^X7X01 — X1X2X4X6X7 + XiX3X 4 X5Xg + X1X2X3X4X6) 

P 6 (x) = l/D 

(X3X4X5X8X9 — X2X3X4X5Xg — X1X3X4X5X8 + X1X2X3X4X5 + X1X2X4X5X6 + 2X1X2X4X5X5X7 
+ X3X5X6X8X9 - X3X5XIX8 - x|x3x|x 6 Xg - X1X3X4X6X8 + X 3 xlx 6 Xsxl - XsX^XgXg 

- X3X4X5X6 + XiX 2 XgX7 — X1X2X3XI) 

P 7 (x) = x 7 + l/D 

I 22,22,22 22 ,32 23 23 22 

t X5X7X8X9 -\- X0X5X7X9 -\- X2X5X7X9 X2X3X5X9 ~T" X5X7X8 X4X5X8 X2X5X7 X4X^X8Xg 

+ x\x$XtX%Xo,+ x\x4X%x\+ x\x 7 x\— x\x\xo,— x\x\x4X%+ x\x4X^X%— x\x\x 7 + 2x\x4X$X§X%Xo> 
— 2X4X5X6X7X8X9 + 2X4X5X6X7X9 + 2X4X5X6X7X8 — 2XIX5X6 — 2X1X2X3X4X5 + X4X 5 X7X8XgJ 

P 8 (x) = x s + l/D 

(24 24 n 3 3 42 24,4 2,22 2,222 

1 X4X6X7 X3X6X8 £X 4X5X5 X4X6Xg X1X4X8 ! X4X8Xg -r X2X4X7X0 t X0X6X8X9 

X2XoX6Xg -j- AX^XqXjXqXq — X -yX 2X 4X7 -\-X-yX 4X7X0 — X4X7X0XQ — ZX\X2X 3 X4Xq-\- AX^X^XQXgXo, 

— 2X4X5X6X7X3X9 + 2X4X5X5X7X8 + 2X2X4X5X6X7X9 — X^X^X^Xg + X2X6X7X9 + X^X^Xs) 

Pg(x) = Xg + l/D 

^42 _ 2 4 „42 ,22 2__222 _2 22,22 2 , q 2 2 

— 2X1X2X3X5X5 + 2x4XgX6X8Xg — 2X4X5X5 + 2XjX 4X5X5X7 'X8 — 2X4X5X5X7X8X9 + 2x4X5XgX7Xg 

24 24,222 222 22 2,22 2 , 4 2 \ 
X3X6X9 X^XqX'j ~y~ Xi X5X7X8 Xi X0X5X8 X5X7X8X0 -r~ X0X5X8XQ -f- XqX-7Xq J , 
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with 



D = D(xi,...,x g ) 



_ 222_222_2 2_222,222_2 2 ,22 

22 ,22 ,22 ,22 222,222 

JuqX'tXqXq -j~ X aXt^X^Xq ~~r XtXcXfXQ -\- XnXi^XyXQ XaX^Xq -\~ XoXaXq 

i 22 ,4 o 222,222 222 

+ X 5 X 6 XrX$ + X 4 X7Xs — 2X1X2X3X4X5X5 — x x x 2 x 3 + x 1 x 2 x 7 — x 4 x 5 x 6 . 

Denote P( k \x) := P^(xi, . . . ,xg) the fc-th iteration of the function P and use the notation 

P^(a,b,c,a,b,c,a + b + c,a + b + c, a + b + c) := (a, b, c, a {k \ b (k \ c (k \ e (fc) , f^, g^A , 
and 

D {k) := D (k) ^b,c,a,b,c,a + b + c,a+b + c,a + b+c) = D(a,b,c, a {k) , b {k) , c (fc) , e (fc) ,f {k \ g {k) ) 

where a^ = a, M°) = 6, c (0 ) = c and e^ = / (0) = ,g (0) =a + b + c. 
Moreover we define, for each id 9 , 

/P 7 (X)4-1 - C fe _! -P4x)lLl -P 5 (x)ltl 

Afc(a;)= -Pi{x)Il_ x P 8 (x)Ik-i - h-i -P 6 (x)I k -i 

\ --P 5 k)^_i -Pe(x)I k -i P 9 {x)h-i - a k -i, 

Theorem 3.5. For each n > 3, the weighted generating function T n (a,b,c) of the spanning trees on the 
Schreier graph £„ of the Hanoi Towers group H^ 3 ' is 



n— 3 

i-n^ (fc) ) a 



T n {a,b lC ) = ^— \[{D^j i " ~ 2 -dct U 2 (P {n ~ 3) (a,b,c,a,b,c,a + b + c,a + b + c,a + b + c))) 

a k=0 



Proof It is clear, from what said above, that T n (a,b,c) ~ C a +b ■ More precisely, the factor a + b 
corresponds to the entry (1, 1) of A„. This entry must be simplified, since we want to compute the 
associated cofactor. If we expand twice the matrix A„, using recursion, and we perform the transpositions 
(17) and (58) for both rows and columns, we get the matrix 

Mn M12 
M 2 i M 22 



I a+b+c 








-d n -2 


— a.I n _2 





-bl°_ 2 








\ 





a + fa + c 


-c-r„-2 





-6/„_2 








— al n _2 










-c/„_2 


a + b + c 








— a In-2 








-bI n -2 




-d n -2 








a + b + c 





-6J„_2 


-al°_ 2 










— al n _2 


-bl n _ 2 








a + b + c 








-d n -2 













— o,I n _2 


-bl„ 2 





a + b + c 








-c/„-2 




-K-2 








-«41-2 





(J 


a + b + c — c T1 _f 


o 










— aI n — 2 








-d n -2 








a + b + c — b n _'_ 







V 





-bIn-2 








-d„-2 








a+b + c — a T1 _ 


2 ) 



Note that each entry is a square matrix of size 3™ 2 and a + b + c is multiplied by I n -2- Hence, the Schur 
complement Formula gives 

(19) det(A„) = dct(A/ 11 )-det(Af 22 -il/ 21 M n 1 Mi 2 ) 



(°^ 3 "- 2 .det 



(D^>) 



-a {1) ll- 2 

~« (1) /°-2 f {1) In-2 b n - 2 







-b {l) ^-2 

-C^I n -2 

g {1) In-2 - a n ~ 2/ 



The fundamental remark is that the matrix in Equation (|19[) has the same shape as A„, since e^°'I n ^i — 
(a + b + c)I n — x- Therefore, we can apply a recursive argument and use the same strategy n — 3 times 
until we get a 9 x 9 matrix, that coincides with A 2 (p( n-3 ) (a, b, c, a, 6, c, a + b + c, a + b + c, a + b + c)) . 

Observe that the entry (1, 1) of A 2 (p(' l ~ 3 )(a, b 7 c,a,b,c,a + b + c,a + b + c,a + b + c)) is still equal to 
a + 6; moreover, all the remaining entries in the first row and column of this matrix are zero. It follows 
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that in the determinant of A2 (p(™~ 3 ) (a, b, c, a, b, c, a + b + c, a + b + c, a + b + c)) a factor (a + 6) occurs, 
that we have to simplify in order to apply Kirchhoff's Theorem. Hence, Equation (|19[) becomes 

det(A„) = (£> (0) ) 3 " 2 -det(A„_i(a, b, c,a, b,c, a + b + c, a + 6 + c, a + b + c-)) 

= (L>(°)) 3 "~ 2 • (D^) 3 "" 3 • det (A„_ 2 (.P (1) (a, 6, c, a, 6, c, a + b + c, a + b + c, a + b + c)j\ 

n — 3 

1 (£)(*0 )3"" 2 . det ( A 2 fp("- 3 ) (a, 6, c, a, 6, c, a + b + c, a + b + c, a + b + c) 

fc=0 

= (a + b)T n (a,b,c). 

a 

4. Some statistics 

In this section, we deal with a statistical analysis about the number of edges, with a fixed label w G 
{a, b, c}, occurring in a random spanning tree of the considered graph. 

It is clear that in the case of the Schreier graphs {S ra }n>i of the Hanoi Towers group, as in the case of 
both directional and Schreier labellings of the Sierpinski graphs {r n } ra >i, all the weights play the same 
role in the labelling of the edges of the graph. 

On the other hand, in the rotational-invariant model, the weights a and b are symmetric, whereas the 
weight c plays a special role in the labelling of the graph. Hence, it is interesting to perform such analysis 
on the Sierpinski graphs {r„}„>i, when the edges are endowed with the rotational-invariant labelling. 

Our techniques are classical: more precisely, logarithmic derivatives of the weighted generating function 
T n (a,b,c) with respect to w give us the mean density of u;-edges in a random spanning tree. We can 
further find the variance and show that the limiting distribution is normal. 

Let w n be the random variable given by the number of edges labelled w in a random spanning tree of 
r„, with w € {a,b,c}. Denote by \x n ,w and <r^ w the mean and the variance of w n , respectively. From 
the remark above it follows that: 

_ 2 _ 2 

fJ-n,a — Mn,b G n,a — °ra,&- 

Proposition 4.1. (1) The means and the variances of the random variables a ni b n and c n are: 

16 ■ 3" + 7 13 • 3" + 1 

2 _ 2 199 ■ 3" + 28 2 _ 34 • 3" - 2 

(2) The random variables a ni b n and c n are asymptotically normal, as n — > 00. 

Proof. Let us prove the assertion for the random variable c n (similar computations can be done for a n ). 
Take the generating function T n (a, b, c) given in Theorem 12.21 and put: 

3"-l 3"+2ra-l 3"- 1 -2n+l 3"~ 1 -1 , . 3" + l 

T n (c):=T n (l,l,c) = 2— 3 J 5 1 (3c+ 2)^^— ( 2c + i)t~. 

We can obtain the mean /u„ iC and the variance a^ c of c„ by studying the derivatives of the function 
log(T„(c)). We get 

/Vc = (log(T„(c)))' | c=1 = ^M c=1 = 13 ^" + 1 . 

T n [c) 30 

Taking once more derivative, one gets 

_ T;'( C )T„(c) - (T r ;(c)) 2 127 -3" + 19 



dog(T n (c)))"|^- (Tn(c))2 

Hence, 
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34 . 3™ _ 2 

r nx = (l0g(T„(c))) | c=1 + Mn , c = — . 
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Next, let C n = — — ^-^ be the normalized random variable; then the moment generating function of C n 
is given by 



3"-3 3" + l 

(133" + l)t 



We get 

E(e to ") = 3 -; ^s Li 5 -: ^'e ! (3«"-2) 1/! f 2 + 3e< 343 "- 2 > 1/2 I I 1 + 2e< 343 "- 2 > 1 

whose limit as n — > 00 is e 2 , showing that the random variable is asymptotically normal. □ 
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